Precise measurements of the 77 and the neutral 
kaon meson masses with the KLOE detector. 
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Abstract 

We present precise measurements of the r] and masses using the processes 
(j) — > ry-y, ry — > 77 and (p KsK^, Ks vr+vr". The mass measurement, 
Mk = 497.583±0.005stat i0.020syst MeV, is in acceptable agreement with the previ- 
ous measurements but is more accurate. We find = 547.874±0.007stat i0.031syst 
MeV. Our value is the most accurate to date and is in agreement with two recent 
measurements based on r/ decays, but is inconsistent, by about lOcr, with a mea- 
surement of comparable precision based on r/ production at threshold. 



1 Introduction 

The ?7-meson mass has changed four times in the past 40 years while the 
measuring accuracy never was better than 0.15 MeV till 2002. In 2005 the 
GEM experiment using the reaction d + p ^ rj ^He at threshold found = 
(547.311 ±0.028stat±0.032syst) MeV [1], while in 2002 the NA48 collaboration 
using the decay r] n^n^n^ found = (547.843 ± 0.030stat ± 0.041,yst) 
MeV [2]. The two results above differ by about eight standard deviations. 
Preliminary KLOE results [3] = (547.822 ± 0.005stat ± COGQ^yst) confirm 
the disagreement. Recently the CLEO-c collaboration found = (547.785 ± 
0.017stat ± 0.057syst) MeV [4] using 4'{2S) rjj/ip decays and combining 
different decay modes of the rj. 

For the K*^ mass there is good agreement between the Novosibirsk (1995) 
and CERN(2002) measurements that have a precision of ~ 30 keV [5,6]. Our 
measurement, similar to that of Novorsibisk with rather increased statistics. 
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is based on the knowledge of the meson mass which is known to 20 ppm 
from the Novosibirsk measurement employing the g-2 depolarizing resonance 
method. The mass is also the basis for the rj mass measurement which 
relies on a precise determination of the collision center of mass energy, W 
in the following. W is determined run by run using e~^e~ —>■ e~^e~ events 
(~ 40, 000 for each run), while the absolute momentum scale is obtained from 
the e'^e~ —>■ (j) KsKl cross section as a function of W. 



2 The KLOE experiment. 

KLOE operates at DA$NE, the 0-factory e~^e~ collider running at a center 
of mass W equal to the (/)-meson mass. Positrons and electrons collide at 
an angle of tt — 0.025 rad. The KLOE detector consists of a 4 m diameter, 
3.2m length drift chamber, DC [7], surrounded by a lead/scintillating- fiber 
sampling calorimeter, EMC [8], both immersed in a axial magnetic field of 
0.52 T with the axis parallel to the bisectrix of the two beam lines. The 
transverse momentum resolution for charged particles is 6p±/p± ~ 0.4%. The 
EMC consists of a barrel and two end caps. Energy deposits in the EMC are 
reconstructed in the calorimeter with energy and time resolutions (Te/E = 

0.057/^^^ (GeV), at = 54 ps / (GeV) in quadrature with 140 ps. The 

centroid of showers is measured with resolution ae = 1 cm/ ^Je (GeV) in the 
coordinate parallel to the fibers and 1 cm in the transverse coordinate. 

For a photon coming from the IP the angular resolution is a ~ lcm/200cm ~ 
5 mrad. Close-by energy deposits are combined into "clusters" . A prompt pho- 
ton is defined as a cluster with \tdu — rdu/c\ < 5at {tdu is the arrival time mea- 
sured at the EMC, is the distance from the e~^e~ interaction point and c 
is the velocity of light) not associated to a charged particle. For this latter, we 
require the distance between the centroid of the cluster and the extrapolation 
of any track reaching the calorimeter to be larger than three times the cluster 
position resolution. 

Only calorimeter signals are used to trigger [9] events for these analyses. We 
require at least two energy deposits above threshold {E > 50 MeV in the barrel 
and E > 150 MeV in the end-cap). The trigger has a large time jitter with 
respect to the event time but is synchronized with the collider radio frequency 
with an accuracy of 50 ps. The time of the bunch crossing producing an event 
is determined off-line during event reconstruction. 

The large cross section for e"'~e~— ~3 fih and for elastic e"^e~ scatter- 
ing, Bhabha scattering, allows KLOE to collect large number of events, some 
450,000 per hour. KLOE takes advantage of these events to maintain a run- 
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ning calibration of time and energy scales of its calorimeter, of the momentum 
and position resolution of the drift chamber, of the machine energy and beams 
crossing angle, and therefore of the center of mass motion, of the mean position 
of the interaction point and of the detector alignment. 



3 Calibration of cm. energy. 

The center-of-mass energy, W, has been measured for each run by fitting 
the e~^e~ invariant-mass distribution for Bhabha events to a Monte Carlo 
generated function, including radiative effects. 

Initial state radiation (ISR), where one or both initial colliding particles ra- 
diate a photon before interacting, affects the e~^e~ collision center-of-mass 
energy W and therefore the final state invariant mass. ISR, which is mostly 
collinear to the beam, is in general not detected. MC Bhabha events were gen- 
erated using the BABAYAGA event generator [10], which accounts for both 
final and initial state radiation. 

An example of this fit is shown in Fig. 1. In a typical run, an integrated lumi- 



entries /.35MeV 




Fig. 1. Fit to reconstructed e^e invariant-mass distribution for Bhabha events, for 
a run with VF=1021.7 MeV 

nosity of 50 nb~^ is collected and W is measured with a statistical accuracy of 
~ 3 keV. The stability of the momentum calibration has been studied measur- 
ing the two-pion invariant mass in Ks — ^ vr+vr" decay. It is found to be stable 
to within 10 keV in the analyzed runs. 
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The center-of-mass energy scale has been cahbrated by obtaining the mass 
from a fit to the cross section measurements for the process e+e~ — >-0 — >■ KsK^. 
The cross section is measured at different values of W around by counting 
the number of Ks tt^ti~ events, correcting for selection efficiency, back- 
ground and the (3%; factor of the KsKl pair, and normalizing to the integrated 
luminosity. Ks tt+tt" events are selected by requiring two tracks with op- 
posite charge that form a vertex within a cylinder of 5 cm radius and 10 cm 
length centered on the interaction point. The invariant mass computed from 
the two pion tracks is required to be within 20 MeV of the nominal neu- 
tral kaon mass. The luminosity is measured by using very-large-angle (> 55*^) 
Bhabha events[ll]. 

The measured cross section is fitted to a theoretical function [12] that depends 
on the (j) parameters, takes into account the effect of ISR, and includes the 
interference with the p{770) and the a;(782) mesons. The mass, total width, 
and peak cross section are the only free parameters of the fit, the p(770) and 
the u;(782) parameters being fixed. The results of the fit to the data are shown 
in Fig. 2. The fitted (p mass is = 1019.329 ± 0.011 MeV, to be compared 
with = 1019.483 ±0.011 ±0.025 MeV, measured by CMD-2 at VEPP-2M 
[13]. The ratio of these two values is used to fix the overall energy scale. The 
correction factor M^^^^ / M^^^^ is 1.00015, corresponding to a shift in the 
value of of ~ 150 keV. 



4 Measurement of the neutral kaon mass. 



The events — »■ KsK^ offer a unique possibility to obtain a precise value 
of the neutral kaon mass. To obtain a crude estimate of the resolution and 
explain the method used we observe that if the (/)-meson is at rest the kaon 
mass can be extracted from the kaon momentum using the relation: 

^k = \^-Pk\ ^— 13 1 

rriK mj^ pk Pk 




Since pk — 110 MeV, measuring it at 1% level, well within the KLOE capabil- 
ity, results in a measurement of the mass better than 0.1%. 50,000 events 
are enough to reach a statistical accuracy of about 1 keV. 

mesons are produced with a momentum along the x axis, = 12.5 MeV at 
DAFNE. From the measured momenta of the two pions from Ks — > ti^ti^ ^ we 
measure the Ks momentum. The Ki momentum is given by px^ = Pep ~Pksi 
where p^ is the average momentum measured with Bhabha events collected 
in the same runs. The center of mass energy of the KsKl pair {WkK) is 
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Fig. 2. Top: Cross section for e^e — KsKl as a function of the center-of-mass 
energy. The soUd line represents the fit to the data. Bottom: Fit residuals. 

related to the kaon mass M^, according to: 



Wkk{Mk) = pMl + 2EksEk, - 2pKs ■ Pk, 



with 

Eks = \/pIs + Ek, = ^pI^ + Ml. 



On the other hand the collision center of mass energy W is computed from 
Bhabha events as described above. 

Corrections due to ISR have to be taken into account when relating W to 
Wkk- The correction function fxiW) has been evaluated using a full detector 
simulation where the radiation from both beams has been implemented, and 
Wkk is reconstructed as in the data. The expression of the radiator function 
has been taken from Ref. [14], including 0{a^) corrections. The correction 
|1 — is very small below the resonance: about 40 keV. It increases up 

to about 100 keV when W is above the mass. In this region radiative return 
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Fig. 3. Mx distribution evaluated in a single run of ~ 16000 events from eq. 2, the 
dotted line is a gaussian fit to the peak. 

begins to be important. The neutral kaon mass is then obtained solving the 
equation: 



The single event mass resolution is about 430 keV. Contributions to the mass 
resolution are: experimental resolution about 370 keV, beam energy spread 
about 220 keV, as measured by KLOE ( in agreement with machine theory ), 
and ISR about 100 keV. The kaon mass distribution for a single run is shown 
in Fig. 3 together with a gaussian fit to the distribution. 

The source of systematic errors considered for this measurement are: 

(1) the momentum calibration; 

(2) the theoretical uncertainty on the radiator function fK^W)] 

(3) the absolute calibration of the beam energy. 

The systematic error due to the momentum miscalibration has been evaluated 
by changing the momentum scale in computing the pion momenta. A momen- 
tum miscalibration, 5p/p , translates to a miscalibration on 5Mk/Mx= 0.06 
5p/p, in agreement with the qualitative calculation made above. The momen- 
tum scale is obtained by using several processes covering a wide momentum 
range 50 to 500 MeV {Kl — > tt^tt^tt^, Kl — vr-^z/, tx^tx'ti^), with a 




(2) 



7 



fractional accuracy below 2 x 10 ^, in agreement with the estimate obtained 
using Bhabhas [15], resulting in a systematic error 5Mk of 6 keV. 

The systematic error coming from theoretical uncertainty of the radiator func- 
tion has been evaluated considering the contribution from higher order terms 
in a. The correction function fxiW) has been evaluated by excluding the 
constant term in the 0{a'^). The corresponding change in fxiW) is 1-3 x 10~^ 
corresponding to a variation on of 7 keV. Further checks have been made 
by using the function given in Ref. [16]: no significant differences were ob- 
served. Additional systematics come from the dependence of the measured 
mass from the W value: we compare the average of the measurements with 
data collected aX W < 1020 MeV with data at > 1021 MeV, where the 
value of fK{W) is more than a factor two larger. The difference between the 
two mass values is Mk{W < 1020) - Mk{W > 1021) = 9 ± 10 keV, consistent 
with zero. 

Other sources of systematics are due to the uncertainties on the W calibra- 
tion, i.e., the statistic and systematic error on 

MCMd-2 obtained 
from our fit. The total contribution from these sources amounts to a mass 
uncertainty of 15 keV. Systematic uncertainties are treated as uncorrelated. 
The result is: 

Mk = 497.583 ± 0.005stat ± 0.020syst MeV. (3) 



5 The 7] mass 

The decay —>■ rj'-f, for </)-meson at rest, is a source of monochromatic ?7-mesons 
of ~ 362.8 MeV momentum, recoiling against a photon of equal momentum. 
Detection of such a photon signals the presence of an ?7-meson. Photons from 
r/— have a fiat spectrum in the range 147 < -E^ < 510 MeV in the labo- 
ratory frame. In the laboratory, the opening angle of the two photons has a 
distribution peaked at its minimum value of 113° . KLOE measures this angle 
with an accuracy of ~0.4°. The value of the minimum angle is a function of the 
rj mass and its measurement determines the mass with a resolution of 2 MeV, 
without energy measurements. In fact we do measure the photon energies. 
The ?7-mass accuracy is however ultimately due to the accurate measurement 
of the photon angles. Together with the stability of the continuously calibrated 
detector and the very large sample of ?7-mesons collected we have been able 
to obtain a very accurate measurement of the ?7-mass [17]. 

Events are selected requiring at least three energy clusters in the barrel calorime- 
ter with polar angle 50° < 6^ < 130°. A kinematic fit imposing energy- 
momentum conservation is performed. The fitted photon energy resolution 
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is vastly improved over the EMC measurement because of the good angular 
resolution. The kinematic fit uses the value of the total energy, the trans- 
verse momentum and the average value of the beam-beam interaction point; 
these values are determined with good precision run by run by analyzing 
e^e~ —>■ e~^e~ elastic scattering events. Fig. 4 shows the of the kinematic 
fit for the data and for Monte Carlo [15] simulated signal events. If more 
than three photons are selected, the combination with the lowest is chosen. 
Events with < 35 are kept for the analysis. 

10^ 



10^ 



10^ 



20 40 60 80 100 120 140 160 180 

Fig. 4. Distribution of the of the kinematic fit: dashed area for the MC simulation 
of — > TT^'j, (/> — > ?77 events; histogram for data. 

Fig. 5 shows the "^7273 '"^7172 Dalitz plot population, with the energies ordered 
as E^^ < Ey^ < Ey^. The m^^^^ ~ m^o, rny^y^ — 'm'^ and rra^^^^ = bands 
are clearly visible. We apply a cut "^^^^^ -|- m^^^^ < 0.73 GeV^, "background- 
rejection cut" in the following, shown by the line in Fig. 5. Events below the 
line are retained for the analysis. The resulting m-y-^y^ distribution, for a data 
subsample, is shown in Fig. 6, top. The 772(7172) distribution in the 542.5 to 
552.5 interval is fitted well with a single gaussian with a = 2.0 MeV as shown 
in Fig. 6,bottom. 

To estimate systematic uncertainties we have studied the effects of the detector 
response and alignment, event selection cuts, kinematic fit and beam energy 
calibration that can influence our measurement. The values of the systematic 
errors are summarized in Table 1. 

First, to check the effect of the e~^e~ interaction point and the alignment of the 
calorimeter relative to the drift chamber, we have selected a high purity sample 
of e~^e~ 7r"'"7r~7 events [18]. The average position of the interaction point, 
determined run by run with e^e~ — e~^e~ events, has been compared with the 
average position of reconstructed tt^tt^ vertex . The difference between the two 
values was computed run by run and the rms of these points (cTvtx) was used to 
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Fig. 5. Population in the "1^273' "^7172 pl^-i^s. The photon energies are ordered as 
E^^ < Ej2 < -£'73- The r] and vr" signal are quite evident. Dashed line see text. 
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Fig. 6. Top: Distribution of the invariant mass 771(7172), for the events selected by the 
cut shown in Fig. 5 Bottom: Distribution of the invariant mass 771.(7172) around the 
value of the 77 mass and the gaussian fit. The result of the fit is ttt,^ = 547. 777 it 0.016 
MeV with xVn-d.f = 168/161, CL=33%. 

evaluate the systematic error introduced in the kinematic fit by varying the IP 
position by ^la^tx- To check for misalignments between the calorimeter and 
the drift chamber, each pion track was extrapolated to the calorimeter and 
compared with the centroid of the cluster. A small correction of 1.1 mm along 
the vertical coordinate, y, and of 2.0 mm along the longitudinal coordinate, z, 
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was applied. The rms (o"dispi) of the mean difference between the coordinates 
of the extrapolated point and the cluster centroid was used to compute the 
systematic uncertainty on the rj mass by shifting the photon point of arrival 
at the EMC by Udispi- 

The energy scale of the calorimeter response and its linearity are checked using 
two different samples of e+e~ 7r"'"7r~7 and e~^e~ — e+e~7 events. The energy 
of the photon, determined from the track momenta and the average value of 
W, was compared with the calorimeter cluster energy. The calorimeter energy 
scale was calibrated to better than 1% and the response is linear to better 
than 2% in the range of interest. The systematic effect on the two-photon 
invariant mass is 4 keV from the energy scale miscalibration and 4 keV from 
the non-linearity. The values above confirm that the mass measurement has 
little sensitivity to the calorimeter energy response. 

It is however important to check the correctness of the position measurement in 
the 24 calorimeter modules of the barrel. We compute the two-photon invariant 
mass for different orientations of the 717273 plane. The rms width of the 77 mass 
distribution was assumed as systematic error: 10 keV and 15 keV respectively 
for variations of the polar and azimuth angle of the normal to the plane. 

Systematic effects from event selection criteria were studied by varying the 
cut (Fig. 4) and the background-rejection cut of Fig. 5. The first has no 
influence (< 1 keV) on the result; the second introduces a systematic error 
of 17 keV determined by translating the solid line in Fig. 5 and rotating as 
shown by the dashed line in the same figure. The mass value is very sensitive 
to the center of mass energy of the 777 system used in the kinematic fit. Due 
to initial state radiation emission (ISR) the available center of mass energy is 
a bit lower than the value computed from the nominal energy of the e~^, e~ 
beams. The effect has been studied with a detailed Monte Carlo simulation 
of the events in the detector, and a shift of ~100 keV was found for the mass 
measurement. Since this correction is relatively large, we have checked the 
MC correction due to ISR emission also for runs taken at different values of 
W. The data were divided in eight energy bins; moreover, two off-peak energy 
bins, centered at W = 1017 and 1022 MeV, were also analyzed in the same 
way as the 0-peak data. Fig. 7, top shows the shift of the mass evaluated by 
MC as a function of W together with the measured shift of the rj mass respect 
to the value obtained a.t W = 1019.6 MeV. Fig. 7, bottom shows the value 
of the mass corrected for the ISR effect. The rms of these points is used as 
systematic error (8 keV). 

The value of the 7r° mass was measured with the same method fitting the low 
mass region of Fig. 6 and the ratio R = mr^/m.^o was also determined. All 
systematic effects discussed above were also evaluated for the mass of the 7r° 
and for the ratio R; the corresponding values are listed in Table 1. 



11 




1025 
MeV 



4 6 
Point number 



Fig. 7. Top. Shift of the value of the r] mass due to ISR as a function of W. Open 
dots: MC, full dots: data (m^(VF)-m^(1019.6)). The line is a fit to the MC. Bottom. 
Value of the rj mass after the correction as a function of W. The band corresponds 
to iblo" around the fitted value. 



Systematic effect 


(keV) 


m^o (keV) 


R (xlO-5) 


Vertex position 


4 


6 


19 


Calorimeter energy scale 


4 


1 


6 


Calorimeter non-linearity 


4 


11 


31 


9 angular uniformity 


10 


44 


120 


(p angular uniformity 


15 


12 


37 


cut 


<1 


4 


13 


Background-rejection cut 


17 


4 


18 


ISR emission 


8 


9 


28 


Total 


27 


49 


136 



Table 1 

Systematic errors evaluated for m^, m^o and the ratio R = m^/m^o. 

Finally, the stability of the results as function of running conditions was 
checked by dividing the data set in eight different periods and determining 
the values of m^, 771,^0 and R for each period. The results are shown in Fig. 8 
and in Table 2 together with their statistical significance. Fit to a common 
value are good. 



The constraint on the total momentum of photons in the kinematic fit is 
very effective in reducing the error on the two-photon invariant mass. The 
absolute scale of W is determined using the CMD2 value as in the K mass 
measurement section. 
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Fig. 8. Value of the r] mass, the vr*' mass and of the ratio measured in the eight run 
periods. The narrow (large) bands correspond to statistical (systematic) error. 





Result of the fit 


xVd.o.f. 


C.L. % 




547,791 ± 7 (keV) 


6.9/7 


44 % 




134,886 ± 12 (keV) 


7.7/7 


36 % 


R 


4.0610 ± 0.0004 


8.9/7 


26 % 



Table 2 

Result of the fit to the values of m^, m^o and R in the eight data periods. 
Rescaling the values of Table 2 for the ratio m'^^^^'^ /m^^'-^^ we obtain 

m^o = (134.906 ± 0.012stat ± 0.048syst) MeV (4) 
= (547.874 ± O.OOTstat ± 0.031syst) MeV (5) 



the vr'^ mass value is in agreement with the world average [19] within lAa. 
As a check of this result, we can use the measured ratio 



R 



m. 



V 



4.0610 ± 0.0004stat ± 0.0014, 



*syst 



(6) 



and the world average value of the 7r° mass, m^o = (134.9766 ± 0.0006) 
MeV [19] to derive = (548.14 ± 0.05stat ± 0.19syst) MeV, consistent with 
the results quoted above although affected by a larger systematic error due 
to a worse cluster position reconstruction of the two photons from ir^ decays 
which have a lower energy. 



6 Conclusions 



Our mass measurement is in acceptable agreement with the previous mea- 
surements shown in Table 3, but more accurate. Averaging [5,6] and our result 
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we obtain M^o = 497.610 ± 0.015 MeV. 



Experiment 


Method 


rriKO (MeV) 


events 


CMD 


[5] 


e+e- ^ KlKs 


497.661 ± 0.033 


3713 


NA48 


[6] 


Kl SvrO 


497.625 ± 0.001 ± 0.031 


665 k 


KLOE 




e'^e^ KlKs 


497.583 ± 0.005 ± 0.020 


35 k 



Table 3 

Recent measurements of the mass. 

Our measurement of the rj (eq. 5) mass is the most accurate result today. It 
is in good agreement with the recent measurements based on rj decays hsted 
in Table 4. Averaging the mass values from [2,4] and our result we obtain 
= 547.851 ± 0.025 MeV, a value different by ~10cr from the average of the 
measurements done studying the production of the rj meson at threshold in 
nuclear reactions. 



Experiment 


Method 


rur, (MeV) 


GEM, MM [1] 
NA48, IM [2] 
CLEO-c, IM [4] 
KLOE, IM 


p d ^He r] 
ri 37r° 
r] 77, 37r^, vr^vr'vr^ 
r/ ^ 77 


547.311 ± 0.028 ± 0.032 
547.843 ± 0.030 ± 0.041 
547.785 ± 0.017 ± 0.057 
547.874 ± 0.007 ± 0.031 



Table 4 

Recent measurement of the r/-meson mass. IM stands for invariant mass of decay 
products, MM for missing mass at production. 
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